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Abstract. Modifications in the production pattern of heavy quark bound states have long been considered 
to provide sensitive signatures for the thermal properties of dense matter created in heavy ion collisions. The 
original concept of Matsui and Satz for quarkonium suppression as signature for deconfinement in heavy 
ion collisions has been challenged recently through lattice studies of spectral functions, which indicate the 
persistence of heavy quark bound states at temperatures well above the transition, as well as through the 
refined analysis of hadronization and recombination models, which take into account the thermal evolution 
of the medium generated in a heavy ion collision. We will review here recent developments on these topics. 
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1 Introduction 

The concept of quarkonium suppression as a probe for the 
thermal properties of hot and dense particle matter or, 
more specifically, for the deconfining nature of the tran- 
sition between hadronic matter at low temperatures and 
a dense partonic medium at high temperatures has orig- 
inally been put forward by Matsui and Satz pQ. In their 
approach they utilized the long-distance infrared proper- 
ties of QCD; in a deconfined medium the heavy quark 
potential gets screened and at sufficiently high tempera- 
ture the screening radius will be smaller than the typical 
size of a quarkonium state. As a consequence the screened 
potential no longer can support the formation of bound 
states. 

Screening strongly depends on a sufficiently high den- 
sity of partons, i.e. color charges that can contribute to the 
screening. In a QED plasma, for instance, the screening ra- 
dius rrj and density n are related through rrj ~ 1/g^/n/T. 
The rapid rise of the parton density at a well defined tran- 
sition temperature, irrespective of this being a phase tran- 
sition temperature or a temperature characterizing a rapid 
crossover, thus intimately connects the occurrence of a 
(phase) transition in QCD with the rapid decrease of the 
screening radius and thus also with the fate of quarkonium 
bound states. 

These qualitative considerations, however, illustrate al- 
ready the difficulty that arises when one wants to use the 
concept of quarkonium suppression as a quantitative tool 
for studies of the QCD transition. On the one hand, it 
is the long distance property of QCD which undoubtedly 
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will prohibit the formation of heavy quark bound states 
at high temperature. On the other hand, as heavy quark 
bound states are small in size also on the typical scale of 
QCD, 1/Aqcd ~ If™, screening has to be strong enough 
to modify also the short distance part of the QCD po- 
tential. Already in the simple context of equilibrium ther- 
modynamics different length scales thus play a role which 
have to be understood quantitatively before quarkonium 
suppression can be established as a unique signature for 
the existence of a phase transition or at least for the char- 
acterization of the properties of the hot and dense high 
temperature phase; a difficulty which is common to many 
signatures for the QCD transition that are currently dis- 
cussed. 

A second ingredient to the concept of quarkonium sup- 
pression in heavy ion collisions which also has been uti- 
lized in pp is that the creation of a heavy quark pair in 
these collisions is a rare event; it is assumed that a gg-pair 
that is created in a dense medium and thus cannot form a 
bound state will separate unhampered and thereafter both 
quarks will no more find another heavy quark as partner 
to form a quarkonium bound state during the stage of 
hadronization. While this was a good approximation for 
the kinematic conditions at the SPS it might not be appro- 
priate anymore at RHIC or LHC. Consequently it recently 
has been suggested that recombination processes |2K| dur- 
ing the cooling of a hot plasma and subsequent pairing of 
heavy quarks initially originating from different creation 
processes (statistical hadronization may lead to a 

different suppression pattern or even to enhanced yields 
of heavy quark bound states. 

A quantitative understanding of the heavy quark po- 
tential at finite temperature, the relevant lengths scales 
that characterize its short and long distance properties 
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and the in-medium properties of bound states combined 
with a thorough analysis of the time evolution of dense 
matter created in heavy ion collisions thus is needed to 
convert experimentally measured heavy quarkonium yields 
into a theoretically well understood tool that can be used 
for the detailed tomography 1 of hot and dense matter. 
In fact, right from the beginning this was understood and 
has led to the discussion of sequential suppression patterns 
for bound states with different quantum numbers its 
momentum and size dependence as well as the role of 
competing suppression mechanisms in ordinary hadronic 
matter jS] . All this has been discussed in detail in a recent 
report on the future heavy ion program at the LHC 

We will concentrate here on a discussion of recent lat- 
tice calculations which aim at a better understanding of 
quarkonium properties at finite temperature. Through the 
detailed study of short, medium and long-distance prop- 
erties of heavy quark free energies on the one hand and 
through calculations of spectral functions from thermal 
quarkonium correlation functions on the other hand much 
better inside into the fate of quarkonium states at high 
temperature has been obtained. In particular, a consis- 
tent picture now seems to emerge from potential model 
calculations, based on calculations of heavy quark free en- 
ergies, and the direct spectral analysis of thermal hadron 
correlation functions. 

We will start in the next section by discussing informa- 
tion on lengths scales relevant for the analysis of quarko- 
nium bound states. Section 3 contains a discussion of the 
calculation of effective potentials from heavy quark free 
energies and their role in potential models. In Section 4 
we present results on quarkonium spectral functions ob- 
tained through a maximum entropy analysis of thermal 
hadron correlation functions. A brief discussion of the con- 
sequences of these results for the analysis of quarkonium 
yields in heavy ion collisions is given in Section 5. Finally 
we give our Conclusions in Section 6. 



1 to use a modern buzzword in heavy ion phenomenology 



2 Short and long distance scales at finite 
temperature 

A quite common qualitative argument for the dissolution 
of heavy quark bound states is that these can no longer 
exist in a hot and dense medium once the screening radius, 
i.e. the inverse of the (Debye) screening mass ro = 1/mjj, 
becomes smaller than the typical size of a quarkonium 
bound state. To go beyond such qualitative arguments it 
is necessary to analyze more quantitatively what small 
and large distances mean in quarkonium physics at a given 
temperature. Some insight into this question can be gained 
by analyzing the short distance properties of excess free 
energies of a static quark anti-quark pair in a thermal 
medium jl()j 

e - -FavO,T) _ Z qq {T, V) 
Z(T,V) 

^\ E ^(TrL(x)Trit (y) ) , (i) 

x : y\\x—y\—r 

or the related singlet free energy, which is obtained from 
static quark anti-quark sources with fixed color orientation 
corresponding to a color singlet state, 

e- F ^T )= l_ J- i(TVL(s)Lt( y )> fl/ . (2) 

x : y\\x—y\—r 

Here L(x denotes the Polyakov loop constructed in lat- 
tice calculations from a product of SU(3) matrices that 
represent the gluon fields |TI]. In Eq. pjjthe subscript, gf, 
indicates that the expectation value has to be evaluated in 
a fixed gauge as the operator used to define F\ is not gauge 
invariant. The corresponding gauge invariant operator, for 
instance in Coulomb gauge, would be a complicated non 
local object |12|. 

The singlet free energy calculated in Coulomb gauge 
is shown in Fig. QJleft) for a selected set of temperatures 
above the transition temperature. The results shown are 
obtained from a calculation in a pure SU(3) gauge theory 
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|13| . Recent results obtained in 2-flavor and 3-flavor 
[T5] QCD are, however, also on a quantitative level quite 
similar. As can be seen, F\(r, T) coincides with the zero 
temperature heavy quark potential, 



4 a 

5 r 



(3) 



for a certain range of quark anti-quark separations and 
then turns over quite rapidly into a constant that de- 
creases with temperature. In fact, these features become 
much more transparent from a calculation of the running 
coupling as a function of distance and temperature [131 
I16| . It can conveniently be defined from a derivative of 
the singlet free energy (T — 0: force) which eliminates 
otherwise uncontrolled constants, 



,(r,T) 



3r 2 dFi(r, T) 
4 dr 



(4) 



This is shown in Fig. fright). At zero temperature the 
Cornell type potential, Eq.[3 yields a qq (r, 0) = a+3or 2 /4, 
with a = 7r/16 arising in string models as correction to 
the linear rising string tension term. This is represented 
in Fig. n by the dashed line at short distances and the fat 
solid line at large distances. Within this parametrization 
of the potential the entire r-dependence of the coupling 
arises from the long distance confinement part of the po- 
tential. Replacing the constant coupling in the Cornell po- 
tential by the perturbatively running QCD coupling leads 
to a characteristic r-dependence also at short distances 
which is indicated by the thin solid line Fig. ^ It gener- 
ally is argued that at T = the perturbative calculation 
of the coupling in the so-called qq-scheme (Eq. is under 
control for r<0.1fm [151117] . 

At short distances the running coupling at finite tem- 
perature coincides with the T — coupling and, of course, 
runs with the dominant scale given by the quark anti- 
quark separation, 



ot m {r,T) 



6tt 



(33 - 2N f )hx(l/rA) 



for rA <C 1 



(5) 



At large distances, however, the coupling defined through 
Eq. 0] drops exponentially, 

a qq (r, T) ~ a{T) exp{-m D (T) r} , for r > 1/T , (6) 

which is in agreement with high temperature perturbation 
theory in a regime where temperature is the dominant 
scale. 

Similar results as those shown in Fig. ^left) for the 
running coupling of the pure gauge theory have been ob- 
tained also in 2 and 3 flavor QCD [Tlll5lll8| . They clearly 
show that (i) the coupling stays temperature independent 
in the short distance perturbative regime for quite a wide 
range of temperatures above T c and (ii) it does so also 
in the non-perturbative regime for temperatures up to 
about 1.5T C ; i.e. the coupling a qq (r, T) continues to show 
remnants of the confining force up to these temperatures. 
Moreover, (iii) the transition to the expected large dis- 
tance behavior, Eq.[S] proceeds quite rapidly at a well de- 
fined distance, r maxi which can be attributed to the point 



r[fm] 
0.8 

0.6 

0.4 

0.2 
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Fig. 2. The scale r me£ j which gives an estimate for the dis- 
tance beyond which the force between a static quark anti-quark 
pair is strongly modified by temperature effects and the Debye 
screening radius, ro = l/mo- Open (closed) symbols corre- 
spond to SU(3) (2-flavor QCD) calculations. The horizontal 
lines give the mean squared charge radii of some charmonium 
and bottomonium states. 



where a qq (r,T) reaches a maximum at fixed T. Due to 
the rapid crossover from short to large distance behavior 
above T c this length scale is, in fact, well be approximated 
by the scale r me <j |14U19j defined as the point at which 
Foo = linv^oo F\ (r, T) equals the zero temperature poten- 
tial V(r). This scale, extracted in a pure gauge theory [TU] 
and 2-flavor QCD ^3], is shown in Fig[2]together with the 
scale set by the inverse Debye mass, ro = I/mo, which is 
extracted from the long-distance exponential screening of 
the singlet free energy. Also included in this figure as hor- 
izontal lines are mean squared charge radii of some char- 
monium (solid) and bottomonium (dashed) states which 
characterize the average separation r entering the effective 
potential in the Schrodinger equation. It is reasonable to 
expect that the temperatures at which these radii equal 
r m ed provides a first estimate for the onset of thermal 
effects in quarkonium states. Of course, as the wave func- 
tions of the various quarkonium states do also reach out 
to larger distances (20j this can only be taken as a first 
indication for the relevant temperatures. The analysis of 
bound states using, for instance, the Schrodinger equation 
will do better in this respect. 



3 Non-Relativistic Bound State Models and 
the Heavy Quark Free Energies 



The influence of a thermal medium on heavy quark bound 
states is commonly discussed in terms of a non-relativistic, 
effective Hamiltonian for the quark anti-quark system G 
I21II22| . Here the influence of the thermal medium is ex- 
pressed in terms of a temperature dependent potential for 
the two-quark system which then can be used to study 
the bound state problem, for instance, by solving the cor- 
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responding non-relativistic Schrodinger equation, 



2m f + — V 2 + V eff {r, T) 

1 TUf 



${ =E((T) $1 



(7) 



where mj denotes the mass of the heavy quark flavor, e.g. 
/ = charm or bottom, i labels different quantum number 
channels and E( is the corresponding binding energy, i.e. 

the bound state mass is then given by M( = 2m / + E( . 
Extensions of this approach taking into account relativis- 
tic corrections have recently been considered in |28j . 

The " effective heavy quark potential" ,V e ff(r,T), mod- 
els the influence of the medium on the qq-system. This is 
then used to solve the bound state problem for a two- 
body system only. The main obstacle in such an approach 
clearly is to obtain a suitable ansatz for V e / / (r, T) . Here 
lattice calculations of the excess free energy of a pair of 
static quark anti-quark sources in a thermal medium, de- 
fined in Eqs. ^and |21 have frequently been used to define 
an effective potential. The often used choice of the color 
averaged free energy as an effective potential, V e /f(r, T) = 
F av (r, T) , may be questionable for two reasons. First of 
all Eq. [7| is meant to solve the bound state problem for a 
color singlet qq-state. F av (r, T), however, gives a weighted 
thermal average over qq-states with arbitrary color ori- 
entation which can be understood in terms of, suitably 
defined, singlet and octet contributions. One thus might 
want to use the free energy for a qq-pair fixed with relative 
color orientation corresponding to a color singlet state, 
V e ff(r,T) = Fx(r,T). Moreover, it might be more suit- 
able to use as an effective potential in Eq. [7| the energy, 
U (r, T), rather then the free energy, F(r, T), of two heavy 
quarks in a thermal bath. The former can be obtained 
from calculations of the free energy using the thermody- 
namic relation, 



U(r,T) 



-T 



2 d F(r, T)/T 
dT 



(8) 



where F and U may correspond to the singlet or color 
averaged free energy and energy, respectively. 



Depending on the ansatz used for the effective poten- 
tial the solution of the Schrodinger equation, Eq. will 
lead to different dissociation temperatures. The qualita- 
tive differences are obvious from Fig. [3] which shows for 
one value of the temperature, T = 1.4 T c , the color av- 
eraged free energy and energy of static quark anti-quark 
sources (left) as well as the corresponding quantities for 
the sources being fixed in a relative singlet orientation 
(right). As the free energy includes an additional entropy 
contribution, F = U — TS, potentials defined in terms of 
the energy are deeper than those defined in terms of free 
energies and thus will lead to higher dissociation temper- 
atures. Likewise, the averaged (free) energies also receive 
contributions from octet configurations, which are repul- 
sive at short distances. This leads to steeper (free) energies 
than in the singlet channel. Dissociation temperatures of 
heavy quark bound states determined from color averaged 
(free) energies are thus generically smaller than the corre- 
sponding values obtained with effective potentials defined 
in terms of singlet (free) energies. Some results for dissoci- 
ation temperatures obtained by using as effective potential 
singlet free energies |23 an d energies j22| > respectively, are 
summarized in Tab.Q] Both calculations are based on lat- 
tice results obtained in the pure SU(3) gauge theory [12]. 
Similar results have also been obtained in |23j . 

For temperatures above the transition temperature the 
free energies and energies that have been calculated now 
also for QCD with dynamical quarks [TlHl3] are quite sim- 
ilar to the pure gauge theory results. The dissociation tem- 
peratures displayed in Tab.^thus give a good description 
of the current status of potential model predictions. As 
may be obvious from Fig. [31 the largest dissociation tem- 
peratures are obtained from XJ\. 

Despite the model dependence entering through the 
choice of V e ff some generic features are common to all 
the potential model calculations. Quite independently of 
the choice of V e /f these calculations suggest that radial 
excitations in the charmonium system (xc, "0') cannot ex- 
ist as bound states in the high temperature phase. In fact, 
the disappearance of % c which will show up experimentally 
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state 


JH> 


Xc 




r 


Xb 




X'b 




El [GeV] 


0.64 


0.20 


0.05 


1.10 


0.67 


0.54 


0.31 


0.20 


T d /T c [from Fi] 
T d /T c [fromf/i] 


1.1 

~ 2.0 


0.74 
~ 1.1 


0.1 - 0.2 
~ 1.1 


2.31 
~ 4.5 


1.13 
~ 2.0 


1.1 
~ 2.0 


0.83 


0.74 



Table 1. Dissociation temperatures in the charmonium and bottomonium system obtained by using pure gauge theory results 
for singlet free energies 1211 (third row) and energies 1221 (last row). 



as a missing contribution to the J/ip yield seems to be a 
good probe for the transition temperature. Calculations 
based on using the singlet energy as an effective poten- 
tial, furthermore, suggest that the ground states (J/ip, ~q c ) 
may still survive in the high temperature phase to tem- 
peratures as high as 2 T c . This, however, is in contrast to 
findings based on effective potentials constructed by using 
free energies rather than energies. This generically leads 
to smaller dissociation temperatures and, indeed, would 
suggest a dissociation of J/ip already at T ~ T c . 

So far we have argued in this section in favor of us- 
ing color singlet free energies (or energies) to discuss the 
fate of quarkonium bound states in the context of poten- 
tial models. Nonetheless, one may also put forward argu- 
ments in favor of an effective potential which incorporates 
some of the color changing processes that take place in a 
thermal medium through the scattering of thermal gluons 
with the static sources. In fact, such processes have been 
discussed as a mechanism for direct quarkonium dissocia- 
tion in a thermal medium |24| . The absorption of thermal 
gluons can change the color orientation of the quark anti- 
quark pair and put it into a color octet state. The analysis 
of color averaged (free) energies thus may have some jus- 
tification also in the context of potential models. 



4 Heavy Quark Bound States from Spectral 
Analysis of Hadron Correlation Functions 

It should have become clear from the discussion in the 
previous section that an analysis of thermal properties of 
heavy quark bound states in terms of potential models and 
probably even the determination of singlet potentials from 
lattice calculations will require additional phenomenolog- 
ical consideration. On the other hand, as pointed out in 
the Introduction, an ab-initio approach to thermal hadron 
properties through lattice calculations exists and is based 
on the calculation of thermal hadron correlation functions, 



G ff (T,r,T) = (J ff (T,r)4(0,0)> 



(9) 



with Jh = q{T 1 r)rnq(T,r) denoting a hadron current 
where the projection onto appropriate quantum numbers, 
H, is controlled by a suitable product of gamma matrices 

The influence of a thermal heat bath on hadron prop- 
erties immediately becomes evident when one compares 
Euclidean time (r) correlation functions of mesons in dif- 
ferent quantum number channels. This is quite apparent 
for light quark mesons where thermal effects are related to 



the weakening of chiral symmetry breaking with increas- 
ing temperature and its restoration at high temperature. 
To illustrate the strong temperature effect on hadron cor- 
relation functions we thus briefly interrupt our discussion 
of heavy quark systems and have a look at the light quark, 
chiral sector of QCD. At low temperature chiral symme- 
try breaking and the breaking of the axial U(l) symme- 
try lead to a splitting of hadronic states, which would be 
degenerate otherwise and would then also have identical 
correlation functions. In Fig. 01 we show pseudo-scalar (ir) 
and scalar (6) meson correlation functions in a gluonic 
heat bath (quenched QCD) at temperatures below (left) 
and above (right) the transition temperature. The spec- 
trum in the pseudo-scalar and scalar channels differs at 
low temperature due to the explicit breaking of the axial 
(7(1) symmetry which consequently leads to quite differ- 
ent correlation functions in Fig. QJleft). In fact, the dif- 
ference will increase with decreasing quark mass values as 
the Goldstone nature of the pseudo-scalar than will lead 
to a flat correlation function. Above T c the almost perfect 
degeneracy of both correlation functions suggests that the 
axial U(l) symmetry is effectively restored. Moreover, the 
significant exponential drop of both correlation functions 
suggests that both modes are " quite heavy" in the decon- 
fined phase of QCD (Fig. Bright)). 

The hadron correlation functions, Gh(t,T), are di- 
rectly related to spectral functions, cth{uj,T), which con- 
tain all the information on thermal modifications of the 
hadron spectrum in the quantum number channels, H, 



G H (r,r,T) 



d 3 p 



(T H (w,p,T) e lp ' 



,cosh(w(r - 1/2T)) 



sinh(w/2T) 



(10) 



From the knowledge of e.g. the vector spectral function 
one obtains, for instance, direct information on thermal 
dilepton rates, 



dW 



5a 2 a v (uj,p,T) 



dujd 3 p 277r 2 c; 2 (e"/ T -l) ' ^ 

Of course, it should be clear that the rates calculated in 
this way do not include any contributions arising from 
the feed-down of other quantum number channels into the 
vector channel |25| . 

Lattice studies of in-medium properties of hadrons have 
greatly advanced in their predictive power through the 
exploitation of the Maximum Entropy Method (MEM) 
|26U27j . This allows the reconstruction of an(uj,p,T) at 
non-zero temperature for light as well as heavy quark 
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Fig. 4. Scalar (8) and pseudo-scalar (it) correlation functions at T = 0.6T C (left) and T = 1.5T C (right). Shown are results 
from calculations with light quarks on a 64 3 x 16 lattice in quenched QCD. 



bound states. We will concentrate in the following again 
on the heavy quark sector and, in particular, will discuss 
in how far the potential model calculations presented in 
the previous section are consistent with studies of charmo- 
nium spectral functions in the high temperature phase of 
QCD 28 , 29 , 3(3| . We will restrict this discussion to the zero 
momentum sector of the charmonium system, <j(uj,T) = 
<jh(uj, 0, T), although also some preliminary results at non- 
zero momentum have been presented recently [311 and also 
some exploratory work for bottomonium has been per- 
formed Eg. 

As indicated in Eq. E3 G H (r,T) ~ / d 3 rG H (T,r,T) 
is related to the spectral functions au(^,T). However, an 
inversion of this integral equation, which would be needed 
to extract <7jy(w,T) unambiguously, is generally not pos- 
sible because lattice calculations only yield information 
on Gh(t,T) at a finite, discrete set of Euclidean time 
steps, r k T = k/N T with k = 0, 1, ...N T - 1, with N T de- 
noting the temporal extent of the lattice. It is, however, 
possible to determine the most probable spectral function 
which describes the calculated data set {Gjy(rfc,T) | k = 
0, N T — 1} and respects known constraints on o~h((*>, T) 
(positivity, asymptotic behavior, ...). This can be achieved 
using a Bayesian data analysis, e.g. the Maximum Entropy 
Method (MEM) gSj- In the context of QCD calculations 
on the lattice this has been introduced in |27|. 

The MEM analysis has its own set of open problems. 
In particular, at finite temperature some ambiguities arise 
through the sensitivity of details of the spectral functions 
to the input spectral function used as default model to 
start the analysis |29II33| . Moreover, also the high energy 
part of the spectral functions is strongly influenced by 
details of the scheme used to discretized the lattice fermion 
action (Wilson doublers 34,35 ). This makes a "point- 
by-point" comparison of spectral functions determined by 
different groups difficult 2 . 

A first feeling for the influence of a thermal heat bath 
on the structure of quarkonium correlation functions and 



2 In any case, one should view the presently available set of 
calculations of spectral functions as exploratory leaving many 
possibilities for improvement. 
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Fig. 5. Thermal correlation functions above T c in axial vector 
(AX: X c ,i) scalar (SC: Xc,o), pseudo-scalar (PS: rj c ) and vector 
(VC: J/ip) channels normalized to reconstructed correlation 
functions which are based on spectral functions calculated at 
T* = 0.75T C (top) and T* = 0.9T C (bottom) 29 . Note the 
different scales on the ordinates of both figures. 



the size of medium modifications of a thermal spectral 
function can be obtained without relying entirely on the 
statistical MEM analysis by comparing directly the nu- 
merically calculated correlation functions at temperature 
T with a correlation function, Grecon(T, T), constructed 
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Fig. 6. Thermal vector spectral functions (J/tp) from MEM analyzes of meson correlation functions calculated in quenched 
QCD on anisotropic (left) and isotropic 29 (right) lattices. The difference in scale on the ordinates of both figures is due 
to different normalizations used for the hadronic currents. 



from the spectral function calculated at a smaller (zero) 
temperature, T* < T, i.e. 



Grecon(r, T) = / duj (T H (y>,T 



cosh(w(r - 1/2T)) 
sinh(w/2T) 

In this way the trivial temperature dependence of the in- 
tegration kernel is taken care of and any remaining dis- 
crepancies between the reconstructed correlation function 
and the actually calculated correlation function at tem- 
perature T can be attributed to changes in the spectral 
function. Some results for charmonium correlation func- 
tions normalized in this way are shown in Fig. [SJ 

The ratio of correlation functions shown in Fig.[S]would 
equal unity, if the spectral functions would not depend 
on temperature. However, as can be seen the correlation 
functions for P-state charmonium (top) , show strong ther- 
mal modifications already at temperatures slightly above 
T c . On the other hand, the 5-state correlators (bottom) 
show only little modifications up to T ~ 1.5 T c . In fact, 
no temperature dependence is visible in the pseudo-scalar 
channel (j) c ) up to T = 1.5 T c and also modifications in 
the vector channel {J/tp) are only of the order of 10% at 
this temperature. 

The reconstruction of charmonium spectral functions 
has been performed using technically different approaches, 
e.g. using point-like [28,29 or smeared hadron sources 
on isotropic or anisotropic |28U;i()j lattice with stan- 
dard Wilson fermion [2H| or improved Wilson fermion |291 
130] actions. The different calculations agree to the extent 
that no significant modification of S-state spectral func- 
tions is observed up to 1.5 T Cl i.e. J/tp and rj c survive 
as narrow bound states with unchanged mass up to this 



temperature. This is shown in the upper part of Fig. 
for J/ip. Current lattice calculations 3 differ, however, on 
the structure of spectral functions for larger temperatures 
(lower part of Fig. EJ). While it is concluded in 28 that 
the J/tp resonance disappears quite abruptly at T ~ 1.9 T c 
the analysis of suggests that the resonance disappears 
gradually; a resonance peak with reduced strength is still 
visible at T — 2.25 T c and finally disappears completely 
at T = 3 T c . 

For the radially excited states (\c) the strong temper- 
ature depends seen already on the level of the correlation 
functions (Fig. [5j also is apparent in the reconstructed 
spectral functions No evidence for bound states is 
found in these channels above the transition temperature. 
Although finite temperature lattice calculations at present 
cannot resolve higher excited states in a given quantum 
number channel, e.g. the tp' in the vector channel, it seems 
likely that also these excitations will get dissolved at T c or 
even earlier. All this is consistent with the revised poten- 
tial model calculations in terms of a singlet free energy. 

To get control over the detailed pattern of dissolution 
of the heavy quark resonances clearly requires more re- 
fined studies. It also should be noted that so far all exist- 
ing lattice studies have been performed in the quenched 
approximation. Although virtual quark loops are not ex- 



3 We stress again that current lattice calculations of spectral 
functions are influenced by lattice cut-off effects which show 
up most strongly at large energies, u>. In fact, only the first, 
low energy peak in the spectral functions shown in Fig. H2 is 
physical and insensitive to changes of the lattice cut-off. The 
other two peaks have been shown to be lattice artifacts arising 
from "Wilson doublers" |34) . 
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SPS 


RHIC 


LHC 


v^GeV) 


18 


200 


5500 


N cB 


0.2 


10 


200 


N ch 


1350 


3250 


16500 



Table 2. Estimate of the number of cc-pairs created in cen- 
tral collisions at the SPS, RHIC and LHC |37| . The estimate 
is based on production cross sections calculated in 39 . Also 
given is the total number of charged particles created in these 
collisions. 



pected to modify the qualitative picture obtained from 
these calculations in the heavy quark sector of QCD one 
clearly has to understand their influence on a quantitative 
level. As discussed in the previous section also the analy- 
sis of free energies in pure gauge theories and QCD with 
dynamical quarks suggests that dynamical quarks do not 
alter the picture gained from quenched calculations signif- 
icantly. 



5 Quarkonium suppression pattern in heavy 
ion collisions 

So far we have discussed possible modifications of quarko- 
nium spectra in equilibrated quark-gluon matter kept at 
fixed temperature. Which impact the basic results of the 
spectral analysis of charmonium, i.e. (i) dissociation of 
excited states close to T c , (ii) persistence of J/tp at least 
up to 1.5T C , will have on observable charmonium yields in 
heavy ion collisions crucially depends on the fate of these 
states during the expansion of the hot and dense medium 
created in these collisions. To this extent it is important to 
understand what will happen to a quark anti-quark pair 
after its initial production in a hard collision during the 
subsequent expansion and cooling of the dense medium 
created around it. If the conditions in the medium are 
such that the initially created pair cannot form a bound 
state and will fly apart, it still may be the case that re- 
combination processes play a significant role during the 
cooling of the hot medium j2BS] ■ If the latter is the case a 
scenario in which the observed charmonium yields may be 
interpreted in terms of a statistical hadronization model 
4 , 5 seems likely to be valid. Contrary to the characteris- 
tic sequential suppression pattern [oTl21| which would re- 
sult from the standard screening scenario for quarkonium 
suppression this would not lead to any threshold effects 
in the energy or temperature dependence of charmonium 
yields. In fact, it may well be that the latter scenario is 
favored for the kinematic conditions at SPS while the sta- 
tistical hadronization picture could become valid at RHIC 
energies and even more so at the LHC. While at the SPS 
at best a single cc-pair is created per collision quite a 
few pairs per collision can be created at RHIC and well 
above 100 pairs could be created in central collisions at the 
LHC. Such an estimate, taken from |3*7ll3*5) . is displayed 
in Tab. 



J/¥/DY 42 . 70 
■ Pb-Pb 
* S-U 
□ p-A 




Fig. 7. Experimental information on charmonium suppression 
at the SPS: Shown is the ratio "measured value" /"expected 
value" for the relative yields B lltl a(J/^p)/a(DY) and 
B'u<r('>p')/(j(DY) as a function of the nuclear absorption 
length L. The "expected value" takes into account the "nor- 
mal" nuclear absorption of charmonium that traverses nuclear 
matter of thickness L. 



As yields in recombination or hadronization models are 
proportional to N^ s , charmonium yields grow much faster 
than the average charged particle multiplicities, 



N 



J/V 



N? £ /N c 



(13) 



At RHIC and LHC the charmonium yields thus could even 
increase rather than decrease irrespective of whether the 
original cc-pair was dissolved by screening or other mech- 
anisms. 

If the formation of bound states is indeed controlled by 
the screening mechanism it may well be that the charac- 
teristic sequential suppression pattern can only be estab- 
lished at the SPS through a refined analysis of the charmo- 
nium system. The currently available data on tp' and J / ip 
yields from the SPS |4THl4*T] . shown in Fig. leave this 
possibility open. In fact, to substantiate the comparison 
of experimental results, which are commonly expressed in 
terms of the nuclear absorption length L 01] , with expec- 
tations based on the screening mechanism it will also be 
important to reach a better understanding of the relation 
between L, the energy density e and the temperature T 
of the dense matter created in a heavy ion collision. The 
relation between L and e is usually established through 
the Bjorken formula, 



dE T /dr/ 
tAt 



(14) 



where the geometrical overlap area At is related to the 
impact parameter and can also be related to the absorp- 
tion length L. Also the transverse energy per unit rapidity 
(dEx /drf) , appearing in Eq. ^] is experimentally accessi- 
ble. Some uncertainty, however, is introduced through the 
choice of the formation time (r) used to characterize the 
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L [fm] e [GeV/fm 3 ] T/T c 

4 ~1 ~1.0 

8 ~3 ~1.16 

10 ~4 ~1.25 



Table 3. Estimate of the relation between nuclear absorption 
length, energy density and temperature in units of the tran- 
sition temperature for Pb-Pb collisions at the SPS assuming 
e c = lGeV/fm 3 . 



initial state of the medium at the time where cc bound 
states might form. This generically is taken to be 1 fm. To 
further relate energy densities to thermal properties of the 
medium it is necessary to know the QCD equation of state 
which is taken from lattice calculations. In particular, the 
energy density at the transition point is presently known 
only with large errors, e c = (0.3 — 1.3)GeV/fm 3 . 

Due to the rapid change of energy density and pres- 
sure in the transition region to the high temperature phase 
|42| small changes in the temperature correspond to large 
changes in the energy density or L. For orientation we list 
an estimate for the relation between these parameters in 
Tab. |31 which is taken from fiures presented in |41) and is 
based on an application of the Bjorken formula, Eq. 1141 
with t = 1 fm. The conversion of this energy density to 
a temperature scale has then been performed using the 
structure of the lattice equation of state by arbitrarily fix- 
ing the energy density at T c to be 1 GeV/fm 3 . This shows 
that the experimentally observed suppression pattern may 
actually correspond to a narrow temperature regime close 
to the transition temperature. 



6 Conclusions 

Quarkonium suppression has been observed in Pb-Pb col- 
lisions at the SPS. The gross features are consistent with 
the original idea of J/tjj suppression arising from screen- 
ing of the heavy quark potential. It is conceivable that the 
suppression pattern observed at the SPS is specific to the 
kinematic conditions in these experiments and that the en- 
hanced creation of cc -pairs at RHIC and LHC may also 
lead to an enhancement of J/tp yields at these colliders. 

Lattice calculations for charmonium spectral functions 
as well as refined potential model calculations suggest that 
the J/tp state can survive up to temperatures well above 
T c whereas excited states dissolve at or close to T c . Given 
the current uncertainties in the conversion of experimen- 
tal estimates for the initial energy density and theoretical 
relations between energy densities and temperatures the 
observed suppression pattern thus is in line with the as- 
sumption that it is caused by the dissolution of ip' and 
Xc states. At least on the theoretical side these uncertain- 
ties can and will be reduced through improvements in the 
determination of the transition temperature and the cor- 
responding energy density. 
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